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Abstract. We investigate the generation of even and odd harmonics using an intense 
laser and a weak second harmonic field. Our theoretical approach is based on solving 
the saddle-point equations within the Strong Field Approximation. The phase of 
the even harmonic oscillation as a function of the delay between fundamental and 
second harmonic field is calculated and its variation with energy is found to be in 
good agreement with recent experimental results. We also find that the relationship 
between this phase variation and the group delay of the attosecond pulses, depends on 
the intensity and wavelength of the fundamental field as well as the ionization potential 
of the atom. 
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1. Introduction 

Attosecond pulses [il |2j are created through the interaction between intense infrared (IR) 
laser fields and atoms or molecules in a process known as high-order harmonic generation 
(HHG). This process is initiated by the creation of an electron wave packet through 
tunneling ionization, followed by acceleration in the laser field and recombination with 
the ion core resulting in extreme ultraviolet (XUV) emission confined to a fraction of 
the laser cycle. The generation process can be controlled by shaping the driving laser 
field, e.g. by coherently adding laser fields of different wavelengths to the fundamental 
laser field [21 IH El E] ■ Most investigations have been carried out by adding the second 
harmonic generated in a doubling crystal [Zl El El [IHl [HI US HSl [H] , with various goals 
ranging from optimization of the conversion efficiency, characterization of the emitted 
attosecond pulses and production of attosecond pulse trains with one pulse per cycle 
of the fundamental field. Several parameters can be varied in these experiments: the 
intensity ratio, the phase difference and the relative polarization direction. 

When the second harmonic field is a weak perturbation to the fundamental, the 
generation of odd harmonics is barely changed, but the induced symmetry breaking leads 
to the appearance of weak even harmonics. High-order harmonic spectra recorded as a 
function of delay between the two fields show that the intensity of the even harmonics is 
modulated and that this modulation has an offset depending on the harmonic order, as 
shown in Figure [l](a) [13j|. These photon spectrograms are at a first glance similar to the 
photoelectron spectrograms of the RABITT (Reconstruction of Attosecond Bursts by 
Interference of Two-photon Transitions) method [H [15] used to characterize attosecond 
pulses. Dudovich et al. [TB] suggested that the two-color HHG spectrograms could be 
used to determine the emission times of the attosecond pulses in situ, thus "probing 
the birth of attosecond pulses". This method was then applied by Doumy et. al. [T7] . 
using laser systems of different wavelengths, and by us [18|, in a direct experimental 
comparison with the RABITT method. 

In this article we investigate HHG by a laser field and its second harmonic by 
solving the Schrodinger equation within the Strong Field Approximation (SFA), using 
the saddle point equation method [TH] . We assume throughout that the second harmonic 
is weak and can be considered as a small perturbation to the fundamental field. We 
calculate the phase variation of the even harmonic oscillation and compare it with the 
harmonic emission times. Our model allows us also to interpret the rapid variation 
of this phase at high energy, observed in several experiments [HI [13], as a change of 
the dominant quasiclassical trajectory from the short to the long one [Figure [l]^b)]. 
The article is organized as follows: In section 2 we review and compare classical and 
quasiclassical electron trajectories for one-color (Sec. 2.1) and two-color HHG (Sec. 
2.2). In section 3 we study the generation of even harmonics close to the cutoff (Sec. 
3.1) and we relate the phase of the even harmonic oscillation to the emission times of 
the attosecond pulses (Sec. 3.2). In section 4 we summarize our results by presenting 
a more general relationship between the phase of the even harmonic oscillation and the 



Probing the birth of attosecond pulses using a two-color field 



3 




Figure 1. (a) Experimental two-color high-order harmonic spectrogram over relative 
delay 13J. The diamonds (o) indicate the experimental delays corresponding to 
maxima of the modulation in the even harmonic orders, (b) Calculated delays 
corresponding to maxima of the even orders for an effective intensity of 1.4 x 10^^ 
W/cm^. Full (dashed) curve are calculated for the short (long) branch labeled 1 (2). 
Below (above) harmonic order 28 the short (long) branch fits the experimental data 
(o)- 



emission times. 



2. Quasiclassical trajectories 

Our method is based on the stationary phase method (also referred to as the saddle 
point method) to solve the Schrodinger equation for an atom exposed to a strong laser 
field using the SFA |20]. In particular, we consider the interaction of an atom with a 
strong IR laser field (frequency u) and a weak second harmonic (frequency 2u). In the 
SFA, the electron will (I) tunnel into the continuum, (II) accelerate in the strong laser 
field and then (III) return to the atom and emit a harmonic photon (frequency Q). The 
phase of the harmonic radiation is related to the quasiclassical action of the electron 

where A{t) is the vector potential of the laser field, Ip is the ionization potential of the 
atom; and p, to, t, m and e are the drift momentum, tunneling time, return time, mass 
and charge of the electron respectively. The high-order harmonic emission will mainly 
originate from the stationary points of 5* — hQt, with respect to all variables [p, t, to] , 
which satisfy the following three equations [20] 

f'dt' eA{t') = {t-to)p, (2) 

ip-eA{to)f _ 
2^ " 
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The electron is thus required to return to the atom at time t [Eq. |2] , to undergo complex 
tunneling at time to [Eq. [s] and to satisfy energy conservation [Eq. |4] . For any realistic 
atom, we have Ip > which implies that the electron must tunnel into the continuum. 
The tunneling process results in damping of the electron trajectories, i.e. complex 
stationary points in the harmonic plateau The electron trajectories beyond the 

cutoff are always strongly damped (also for Ip = 0), because they are always classically 
forbidden, corresponding to large imaginary components of the stationary points. The 
physical emission time of a given harmonic is given by the real part of the complex 
emission time [151 l22l 123] . 



2.1. One-color case 

We first consider the one-color case with a vector potential, A(t) = Aism{u)t), where 
u is the angular frequency of the fundamental laser light. The complex stationary 
points [p<"\t("),tj''^] are calculated as a function of the high-order harmonic photon 
angular frequency, Q. The index n is used to separate different sets of solutions, where 
n = 1 corresponds to the short branch; and n = 2 corresponds to the long branch. A 
direct comparison between the stationary points for Ip = 15.76 eV (Ar) and Jp = eV 
(referred to as the classical case) is shown in Figure [2] The detailed behavior of these 
stationary points will prove to be important not only for the generation of attosecond 
pulses from one-color HHG, but also in the quantitative analysis of two-color HHG. The 
classical case corresponds to solving Eqs. [2]|4] for Ip = 0. The stationary points are first 
real, describing classical trajectories in the continuum. Using classical trajectories to 
explain the HHG process leads to the approximate cutoff law: hQ^ax ~ 3.2Up -\- Ip 
[2^ |25] . Harmonics above this cutoff can not be generated because there are no 
electrons returning with suffiecient kinetic energy. Quantum mechanically, we see that 
the stationary points quickly develop large imaginary components above the cutoff, 
which result in an exponential damping of the harmonic yield. Solving Eqs. [2]j4] for 
Ip 7^ 0, the exact cutoff law is found to be fiVtmax = 1-3/p -|- 3.2[/p [20] • In Figure |2] 
we have, therefore, shifted the classical solutions by 1.3/p so that the cutoff harmonic 
of the classical case coincides with the quantum mechanical case of argon for a more 
meaningful comparison. The similarities and differences between the classical case (red 
line) and the quantum case (blue line) are clearly observed: The real parts of the drift 
momenta [Fig|2](a)] are almost the same in the two cases, apart from a small systematic 
shift. The real parts of the return times [Figjlj^c)] overlap almost perfectly in the 
main part of the plateau for both the short (line) and long (dashed) branches. This 
demonstrates the usefulness of the classical model for understanding the intrinsic chirp 
of the returning electron wave packet and the resulting chirp of the attosecond light 
pulses. The real parts of the tunneling times [Figj2|^e)] are, however, quite different. 
In the lower part of the plateau there is a discrepancy of a factor of two for the short 
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Figure 2. Real and imaginary parts of the stationary points, [p,t,to]. Ionization 
potentials: argon ij/^^'' = 15.76 eV (thick blue) and classical /p°'' = eV (thin red). 
The classical case is shifted by 1.3/j, to match the quantum mechanical cutoff. The 
short branch is a line; while the long branch is a dashed line. The intensity is 2 x 10^^ 
W/cm^ and the wavelength is 800 nm. 



branch. The short branch trajectories tunnel at earher times compared to the classical 
case due to the ionization potential. This implies longer times in the continuum and 
an increased amount of phase acquired, as the path of integration extends in Equation 
[ij For the long branch the trend is opposite and smaller. The imaginary parts of the 
stationary points [Fig|2]^b,d,f)] lead to exponential damping of the electron amplitudes. 
This is especially clear for the imaginary part of the tunneling time [Fig|2](f)] where the 
complex part is zero in the classical case where no tunneling is required; while, in the 
complex case, the short branch will suffer more damping than the long branch due to a 
lower instantaneous electric field strength at the time of tunneling. 

The high-order harmonic dipole from one half period of the laser field can be 
approximated as a sum of stationary contributions 

r rT/2 ft 

= j d^V j dt j dtoAexp[iS{p,t,to)/h-iQt]^ ■ (5) 



n=l,2 

These discrete contributions correspond to different sets of stationary phase solutions 

1/2 

X 



'n = A^"^ I I exp ks^-^h - zfit^")] , (6) 
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where S^""^ = 5(p<"\ t^"), and where 

Mij = — [S{p,t,to)-hQt], (7) 

is the Hessian matrix of second-order derivatives of the Legendre transformed action with 
i and j being any of our variables [p, t, to] . The Hessian matrix is complex symmetric 
and we have assumed that it can be diagonalized so that the one-dimensional stationary 
phase approximation can be applied to all five integrals independently. The prefactor 
A = d{p — eA{t))* E(to)d{p — eA(to)), describes the dipole transitions to the continuum 
at time t^, and the subsequent recombination to the ground state at time t. In the 
following, we will assume that A is slowly varying and that it, therefore, does not affect 
the stationary points. A more exact analysis could, however, also include the effect of the 
energy dependent atomic scattering phases, rji, in the stationary phase equations (Eqs.[2]- 
|4| (from the recombination matrix element d(p — eA(t))* oc exp[—irji]). This would lead 
to a small change in the stationary points, which would be especially interesting to 
study in argon due to an unusually strong variation of the scattering phase in the 
present energy region |26j . 

Assuming the process to be periodic, the total dipole response, Xq, is found by 
summing the stationary points from one whole period of the fundamental laser field, 
< t < T, i.e. two adjacent half periods. There are thus two contributions from each 
branch, 

(t + r/2, to + T/2) = -4") (t , to) exp [-^^]^/2] , (8) 

separated by a half period, T/2. The overall minus sign comes from the change of sign 
of E{t); while the phase factor originates from the Fourier component in the Legendre 
transformed action. The total dipole response becomes 

x„ « E f { 1 - exp i-inm } = 2 1: 4"' X 1 1; "^'^ . (9) 

„ n [ 0; "/i^ IS even 

where the two contributions add constructively for odd orders; and destructively cancel 
for even orders. 

2.2. Two-color case 

We now consider a two-color laser field composed of a fundamental laser field and a 
weak second harmonic with the same polarization, A{t) = Ai sm{ujt) -\- XA2 sin[2ut -\- (p) , 
where A is a small perturbation parameter. It is possible to solve the two-color high- 
order harmonic emission using the stationary phase equations [Eqs. |2||4] directly, but 
this requires evaluation of the system at all values of 0. We will follow a different route 
where the second harmonic is treated as a perturbation and only the stationary points 
of one-color HHG need to be calculated. The two-color action is expanded in A as 

s^r,4i^iM]i^r\,f',,i^iMmMM, (10) 



to \ 2m / Jto m 
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The first term in Eq. [TO] corresponds to tlie action in the one-color case, Si, and the 
second term is the correction term due to the interaction with the second harmonic field, 
a = The correction term can be calculated as 



e 

a = — 
m 



^-^ cos(2ut' + 0) r r / sin(3wt' + (f)) sin(ut' + 
7: \- eAiA2 



,(11) 

to 



2uj \ Quj 2uj 

where a depends on the ionization potential through the stationary points. Note that the 
derivation of Eq. 11 using the quantum mechanical stationary points differs from those 
presented previously [HI [18] , because it includes the effect of the ionization potential 
within the SPA. 

The two-color high-order harmonic dipole can be calculated by using Eq. ^ and by 
adding an additional slow factor due to the second harmonic. In analogy with Eq. 
the dipole from one period of the fundamental field contains one discrete contribution 
for each branch and half period. 



1 exp 



2E 



(n) 




— exp 



(n) 



inr 



2h _ 

fl/u is odd 
Q/u is even 



(12) 



where Xq^ is the half period contribution in the one-color field, and where the following 
property 



^(") (t + T/2, to + T/2) = -a^"^ (t, to 



(13) 



relating the two-color phase between adjacent half periods for a given branch n is 
used. The intensity of the even order harmonic emission from a single atom can be 
approximated by 

2 



In oc 



Xc 



(14) 



which is valid for \a'''^^ \ <^ n. We define the intensity from a specific branch, n, as 

4"H0)oc|a(")(0)f. (15) 

This (artificial) separation of the branches can be realized in macroscopic medium either 
by phase matching in a long gas cell [T3| or by spatial separation in the far field [HI [TS] . 



3. Subcycle delay dependence of even harmonics 

One fascinating aspect of two-color HHG is that it depends on the subcycle delay (or 
relative phase 0) between the laser fields. Dudovich et al. [H] have proposed to use the 
(/)- dependence of even-order harmonics in order to estimate the emission times of the 
attosecond pulses in situ. We will refer to the relative phases that maximize the intensity 
of the even harmonics as the in situ phases: (j)Q\Vt) and (/)o^^(f2), for the short and long 
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branch respectively. The in situ phases are plotted in Figure |3] (a) and (d) as a function 
of harmonic order for argon at two different IR intensities. The corresponding emission 
times of the high-order harmonics, Re{t^"'\Q)} = t^^(fi), are shown in (b) and (e). The 



40 



35 



30 



03 



25 



(a> 


1 


1 


1 


\ 




\ 






! \ 












» 1 



(b) 






I, 










/ 









0.5 

Delay [T/2] 



(c) 


























Emission time [T] 



0.5 

Delay [T/2] 



(f) 


1 














V 

















Figure 3. Quantum mechanical in situ phases, (/'o"\ between the fundamental 
and weak second harmonic field that maximize the even harmonic emission in argon 
(4^'') = 15.75 eV) for intensities of (a) 1.5 x 10^^ W/cm^ and (d) / = 4 x lO" 
W/cm^. The in situ phases are presented as delays, Tq"' — — (/)q"''/2w in units of T/2. 
The corresponding high-order harmonic emission times are shown in (b) and (e), in 
units of T. Classical in situ phases (/p°'' = 0) for intensities of (c) 1.5 x 10^** W/cm^ and 
(f ) / = 4 X 10^'' W/cm^. The short (long) branch is a full (dashed) curve. T = 2tt/u} 
is the period of the fundamental laser field. 



intensity in (a,b) is realistic for typical laser pulses of 30 fs in argon; while the intensity in 
(d,e) is greater than the saturation intensity and, therefore, not experimentally feasible 
for multicycle laser pulses. The aim of our analysis is to distinguish between what can 
be measured experimentally, 0o"'''(f2), and the desired emission times, t^^\Q). The in 
situ phases from the classical model (/p = 0) are plotted in Figure |3] (c) and (f) for 
comparison. 



3.1. Behavior close to cutoff 

The similarity between the in situ phases and the emission times is striking, especially 
for the high intensity [Fig. [3](d,e,f)] where the short and long branch merge in the cutoff 
at harmonic 63. At lower intensity [Fig. [3] (a)] the in situ phases do not merge in the 
cutoff. Intuitively, one might think that the short and long branch should merge in 
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the cutoff, but this is not necessary in the quantum mechanical case. The stationary 
points of the short and long branch do not merge on the imaginary axis in the cutoff 
[Fig.|2](b,d,f)]. The different behavior of the in situ phases for the short and long branch 
in the cutoff is, hence, an amplitude effect rather than a pure phase effect. This is verified 



by inserting only the real part of the stationary points into Eq. (15) which does indeed 
yield coincidental cutoff behavior of both branches, as expected for pure phase effects. 
Furthermore, it is the long branch that remains physical beyond the cutoff; while the 
strange behavior of the short branch arises from a set of stationary points that become 
unphysical beyond the cutoff [27j. In the classical model [Fig. [s] (c,f)], the in situ phases 
always merge in the cutoff since there are no amplitude effects (damping) in the plateau 
[Fig.|2|(b,d,f)]. 

It is also worth to note that the in situ phases of the two branches intersect at 
lower harmonic orders than the cutoff. This "intra-plateau crossing" is marked with a 
cross ( X ) in Figure [3] (a,d) and it should not be confused as the position of the cutoff. 
At high intensities it is easy to distinguish between the intra-plateau crossing and the 
cutoff; while at low intensities they may be separated by a few harmonic orders only. 
In Figure [T] we present an experimental result where 0o is linear from harmonic 22 to 
harmonic 28. At higher harmonic orders, a dramatic bend is observed [T8] . Using 
our quantum mechanical model we identify the lower orders as part of the short branch, 
n = 1, (full curve); while the higher-orders are identified as the long branch, n = 2, 
(dashed curve). We stress that the dominance of the long branch close to the cutoff 
is a new result which appears already at the single atom level when the short branch 
becomes unphysical. This new finding illustrates the usefulness of quantitative probing 
of HHG using a perturbative field. This effect would be very difficult to observe with 
the RABITT method because the corresponding emission times always merge in the 
cutoff, see Fig. |3] (b,e). The classical model fails to reproduce the bend [Fig. |3] (c)] 
because the corresponding intra-plateau crossing occurs at much lower harmonic orders 
(not shown). Having discussed the details of the behavior close to the cutoff, we now 
turn our attention to the central part of the harmonic plateau. 

3.2. Ratio of in situ phases and emission times 

In the following, we consider the first spectral derivative of the in situ phase, d(l)Q^'^ /dfl, 
which we wish to compare with the group delay dispersion (GDD), dtP/dn. Both 
quantities are evaluated in the central 50% of the harmonic plateau, i.e. in a region 
set by 1.3/p + 3.2Up ± 0.25 x 3.2Up, where both and t^^^ are linear to a very 

good approximation. We avoid fast and nonlinear variations both close to the ionization 
potential {hfl = Ip) and close to the cutoff {hfl = 1. 3/^ + 3. 2f/p) using this central region. 
We define the ratio between the two quantities as 

= -uj^ = -uj^^ (16) 
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The GDD can in principle be obtained from the 0-dependence of the even-order 



harmonics as dt^^^ /duj 



(n) 



/du/u. This relation would be very useful if 7 



(n) 



was a constant (or at least a constant for each branch n). Unfortunately we show that 
7^"^ depends on the ionization potential of the atom, as well as the laser intensity and 
wavelength. 

3.2.1. Role of ionization potential In this subsection we study how 7*-"^ varies with 
laser intensity for a given laser wavelength (800 nm) and different atomic species (helium, 
argon, sodium and the classical case). In our model of 7^'^-', atom specific properties 
enter only through the value of the ionization potential, Ip. In Figure |4| we plot 7*^"-* 
for the short branch (a) and the long branch (b). In the classical case {Ip = 0) we find 
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Figure 4. Ratios for different ionization potentials as a function of fundamental 
intensity for (a) the short branch and (b) the long branch. Symbols are o helium 
4""^ = 24.58 eV, □ argon 4^"^^ = 15.76 eV, V sodium 4^"' = 5.14 eV and o 
classical Ip^^ = eV. 

constant ratios: 7^^-' ^1.1 for the short branch and 7*-^'' ~ 0.84 for the long branch; 
while in all realistic cases, where Ip > 0, there is a significant change in 7^"^ as a function 
of the fundamental intensity. In the case of argon there is an especially strong variation 
of the ratio: from 1.6 to 0.5 for the short branch over the given intensity range. This 
strong intensity dependence is problematic for the experimental determination of two- 
color HHG since the exact effective intensity is often unknown in experiments. The 
variation of the ratios is smaller for the long branch, because the stationary points are 
more similar to the corresponding classical case, as seen in Figure [2j The variation 
of the ratios is smaller for sodium than for helium because of the smaller ionization 
potential. Sodium is chosen as an example to show that there is a significant difference 
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between the classical and quantum mechanical cases even for a relatively small ionization 
potential. To use the two-color HHG approach for characterization of attosecond pulses 
would require an accurate determination of the laser intensity, as well as a quantum 
mechanical calculation of 7'-"^''. 

Next we comment on our previous work [IB] where we compared the in situ method 
and the RABITT method for GDD of attosecond pulses. An explanation for the good 
agreement we found between the two methods may come from the crossing of 7*-^^ for 
argon occurring at 1.8 x 10^^ W/cm^ with the classical limit, see Fig. |4](a). It would be 
interesting to see more experimental results, carried out on different atoms. 

3.2.2. Role of laser wavelength Finally, we study the dependence of the ratios, 7^"^ 
with laser wavelength. The ratios are calculated for 800 nm, 1.3 /im and 2 /im, 
corresponding to a Titanium-Sapphire laser system and two mid IR laser sources. We 
use the ionization potential of argon. Similar to Fig. |4| all ratios look different when 
they are plotted as a function of laser intensity, see Figure |5] (a) for the short branch 
and (b) for the long branch. The ratios of the short branch show larger variations than 

2 
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1.5 

1 



0.5 

1 1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 3 

Intensity [W/cm^] xio^* 

Figure 5. Ratios for different laser wavelengths as a function of fundamental intensity 
for (a) the short branch and (b) the long branch. Symbols are: o A = 800 nm, □ A = 1.3 

(Ar) 

/im, A A = 2 /im are for argon Ip = 15.6 eV; while symbol o is the classical case 
Ip°^ — eV, which is independent of A. 

those of the long branch. Furthermore, it is seen that the longer wavelengths lead to 
less variation of the ratios over intensity. The recent experiments of Doumy et. al. [T7] 
were carried out at these wavelengths. Applying our method improves the scaling law 
for the GDD from A"'^'^'' to A~^ °^, which is closer to the expected scaling of the 
harmonic chirp times intensity. We stress that in order to use the calculated ratios to 
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determine the GDD, the experimental intensity must first be measured as accurately as 
possible independently of the two-color HHG scheme. 



4. Discussion and conclusions 



In our quantum mechanical derivation of two-color HHG we find that both amplitude 
and phase effects are important. The ratio of the ponderomotive energy and the 
ionization potential, Up/Ip, serves as a measure on how "classical" or how "quantum" 
the electron trajectories are. Choosing this ratio as our x-axis for the data in Fig. (|4])-([5|, 
we find that all individual ratios 7*^"\ follow a universal curve, as shown in Figure [6j 
Photoelectron emission ranging from the photon picture to the tunneling picture is 
described in the theory of Keldysh |28] where the ratio ^Jlp/2Up <^ 1 (corresponding 
to Up/Ip 3> 1) implies efficient tunnel ionization. In this limit we find that 7^"-* slowly 
converges toward the classical limit, Ip = 0. One should look at Fig. [6] with some 
caution since the SFA is derived for the long wavelength limit requiring: Ip ^ hu and 
Up/Ip > 1/2. It is clear that the behavior of the HHG process changes dramatically 
around Up/Ip ~ 1, i.e. when the kinetic energy of the electron is close to the potential 
energy of the atom. 




Figure 6. Universal ratio between the group delay dispersion and the derivative of 
the in situ phases as a function of the ratio of ponderomotive energy and ionization 
potential for (a) the short branch and (b) the long branch. The classical limit, Ip^^ = 0, 
is reached very slowly see gray arrows. 



In conclusion, we have studied the HHG process perturbed by a weak second 
harmonic field within the SFA. We find that the dependence of the even harmonics on 
the subcycle delay between the two fields, can not be understood using classical theory. 
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Our calculations show good agreement with experimental results [T3l [18] showing the 
change of behavior at high energy, explained as a change of dominant quasiclassical 
branch of trajectories. We stress that there is an intra-plateau crossing between the 
short and the long branch which does not coincide with the true cutoff. Furthermore, 
we calculate the ratio between the GDD of the attosecond pulses and the phase variation 
of the even-order harmonics, 7*^"\ as a function of intensity, wavelength and ionization 
potential. The analysis method called in situ probing of the birth of attosecond pulses 
[16] must be improved by considering the influence of the atomic properties and laser 
parameters, before it can be applied for quantitative experimental studies. Using the 
classical analysis will only lead to a qualitative prediction for the GDD of the attosecond 
pulses with the correct sign. It will be interesting to compare the results of our quantum 
mechanical approach to probing the birth of attosecond pulses using a two-color field 
with more refined calculations. 
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